PERIODIC LATTICE WITH DEFECTS 



ANTON KUTSENKO 



Abstract. The discrete periodic lattice of masses and springs with line and point defects is 
considered. The dispersion equations for propagative, guided and localised waves are obtained. 
" The detailed analysis of example with three masses is provided. 
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There are many papers devoted to the wave propagation through the discrete periodic lattice 
with different types of defects. One of the early works [M] provides explicit solution of the 
wave equation for the periodic lattice with two different masses. The structure with an infinite 
line defect embedded in an infinite square lattice has been considered in the paper |OAj . In 
Ph| this case the dispersion relations for localised modes can be computed in explicit form. In the 

paper [MS] the authors examined several classes of continuous and discrete models with various 
defect configurations. In the recent paper jCNJMM] the homogeneous lattice perturbed by 
the finite line of masses is considered with respect to analysis of eigenmodes. Also note the 
large number of papers devoted to periodic structures without localised defects but with the 
boundaries, see references in |KKlj . 

The motivation for the present paper is to combine perturbation of the periodic lattice by 
a periodic line (like in |OAj ) with that by a finite defect (like in |CNJMM] ) and to obtain a 
new dispersion equation for such configurations. The developed method can be applied for 
the different types of lattices with arbitrary horizontal and vertical periods. Note that this 
method is based on the modified Green function approach which differs from the method of 
|M1] . Further we plan to adapt the so-called monodromy matrix approach (see e.g. [KS]) 
for obtaining dispersion equations for localised modes. In particular the monodromy ma- 
trix approach should allows us to obtain localised modes in propagative and guided spectral 
intervals. 

In the example below, we consider the homogeneous lattice with mass of particles M = 1 
perturbed by the line defect with mass of particles M and with one single mass M. The de- 
tailed analysis of the band-gap structure of propagative and guided spectra and the necessary 
and sufficient conditions of existing of localised modes are obtained. We note an interesting 
phenomenon: if the mass M increases indefinitely then the upper limit of mass M for which 
the localised mode exists is exactly j — j- (somewhat unusual appearance of n). 

The work is organized as follows. Section [5] contains our main results, namely, the dispersion 
equations for three types of defects: 1) periodic lattice without defects 2) periodic lattice with 
periodic strip 3) periodic lattice with periodic strip and with finite defects. In Section [3] 
the derivations of dispersion equations from the Section [2] are provided. In Section |4] the 
propagative, guided and localised spectra for the homogeneous lattice with line and single 
defect are analised. 
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2. Dispersion equations 




Figure 1. Periodic lattice with defects. 



2.1. Periodic lattice. Consider the 2D periodic lattice (see Fig [TJ without defects, i.e. M n = 
M n = M n . The wave equation takes the following form 

u « - 4n " = -^ 2 M n u n , (2.1) 

n'~n 

where u is the frequency, u n is antiplane displacement and n = (ni,n2) is a position of the 
mass M n on the square lattice (number of the nod). We suppose that all quantities in (12. ip 
are normalised (dimensionless). The notation n ~ n' means that there exists a link between 
nodes with numbers n and n'. Suppose that the function M n is periodic 

M n+Niei+N2B2 = M n , VnGZ 2 (2.2) 

with basis vectors e x = (1,0), e 2 = (0,1) and periods Ni,N 2 ^ 1. After applying Fourier 
transformation (see (j3.4p - fl3.6D ) the equation (12. ip takes the following form 

Lu = -u 2 Mu (2.3) 

with finite matrices defined on the unit cell M = [l..Ni] x [l..jV 2 ]: 

M = (M n 5 nn /) nXeAr , L = L -4i, (2.4) 

Lo = (^n~n')n,n'eA^5 I = t^nn' )n,n'6A^; (2-5) 

where <5 nn / is a Kronecker symbol and 

5 n ~n' = 5n~n'+ (2.6) 
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^— iki r i „iki c , 

e Cn-Afiei~n' + & On+Ar 1 ei~n' + 

e On-A^'-n' + e O n +7V 2 e 2 ~n' 

with 5 n ^ n ' = 1 iff there exists a link between nodes n and n' and 5 n ~n> = otherwise. 

The dispersion equation which determines the Floquet branches cu p (k) (dispersion curves for 
propagative spectrum) can be written in the form 

det Lp(w, k) = 0, where L p = L + cu 2 M. (2.7) 

2.2. Periodic lattice with periodic strip. Here the periodic strip with masses M n = 
M n + M^ 1 ' (see Fig. [1]) is added to our lattice. In this case along with the propagative 
spectrum u p there appears the guided spectrum u g . The guided spectrum corresponds to the 
waves which are bounded (quasiperiodic) along the added strip and decay along the direction 
which is perpendicular to the strip. The guided dispersion curves depend on Floquet parameter 
ki and can be determined from the dispersion equation (see the derivation in Section | 



detL g (w, fci) = where L g = i + w 2 (Lp l ) 2 Mj (2.8) 

and 

Mi = (M^5 nn ,) n , n/eM , (•), = -!- / -dkj. (2.9) 
Note that the equation ( 12. 8 p is not valid in the bands of propagative spectrum, i.e. in the sets 

I p (A;i)=Wp(A;i,[-7r,7r]) ( 2 - 10 ) 

which correspond to the projection of the propagative dispersion curves w p (k) on the plane 
(u,ki). We need this restriction because the inverse of L p (see (12.81) ) does not exist in the 
intervals I p (ki). 

2.3. Periodic lattice with periodic strip and with localised inclusions. Here the lo- 
calised inclusions with masses M n = M n + ML 2) (see Fig. HJ are added to our periodic lattice 
with the strip. In this case along with propagative u p and guided u g spectra the localised 
spectrum <jj\ oc appears. The localised spectrum corresponds to the waves which decay in any 
direction of our lattice. The localised spectrum u\ oc can be determined from the following 
dispersion equation (see the derivation in Section 13. 2p 

AocM = det(i + oo 2 (L g 1 (L p 1 ) 2 ) 1 M 2 ) = 0, (2.11) 

where 

M 2 = (M( 2 )<W)„X6Af- (2-12) 

Note that the equation (12.111) is not valid in the bands of propagative and guided spectra, i.e. 
in the sets 

Jp = / P ([-7r,7r]), J g = a; g ([-7r,7r]) (2.13) 

which correspond to the projection of propagative w p (k) and guided u) g (ki) dispersion curves 
on the axis u. We need this restriction because the inverse of L p and L g (see (12. lip ) do not 
exist in the intervals I p and J g . 
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3. The derivation of dispersion equations 

3.1. Periodic lattice. The set of linear equations (12. ip can be represented in the form of 
infinite matrices 

Lu = -w 2 Mu, where u = (w n ) neZ 2 (3.1) 

and 

M=(M n <W) n ,z,' 6 z», L = L„-4I, (3.2) 

Lo = (^n~n')n,n'eZ 2 ? I = (<W )n,n'eZ 2 • (3-3) 

The infinite system (13.1 ft can be rewritten in the Fourier space as a finite system. For this we 
introduce the operator 

^:f(Z 2 )^([-7r,7r] 2 ), 7Xu)=u(k), (3.4) 

where 

u(k) = (u„(k)) neA /-, tt n (k) = ^2 Un+r 1 N 1 e 1 +r 2 N 2 e2 ex P 

(ir-k). (3.5) 

rez 2 

Under the action of the operator J 7 any sequence u from ^ 2 (Z 2 ) becomes the vector-function 
u(k) with N1N2 components, any component is a function from L 2 ([— tt, tt] 2 ). Then the wave 
equation ( 13. ip can be rewritten in the form 

JXJ^u = -u^JTVLF^u (3.6) 

which coincide with (12. 3p because J r LJ r_1 = L and J r MJ r_1 = M. 

3.2. Periodic lattice with periodic strip. In this case the wave equation (I2.3P becomes 

Lu = -u 2 (Mu + Mj (u) 2 ) (3.7) 
with Mi defined in ( 12. 9p . From the equation ( 13. 7p we obtain that 

u=-w 2 Lp 1 M a (u) 2 (3.8) 
with L p defined in (12. 7p . Integrating with respect to k 2 we get 

(u) 2 = -w 2 (L p 1 ) 2 M 1 (u) 2 (3.9) 

or (see definition of L g in ( 12. 8p ) 

L g (u> 2 = 0. (3.10) 

From ( I3.10P we obtain the condition of existence of guided waves ( 12. 8p . 
For our purpose we need the following fact. Consider the equation 

Lu + w 2 Mu + w 2 Mi(u} 2 = f (3.11) 

with some vector-function f from Lj^([— tt, 7t] 2 ). Then the solution of this equation takes the 
following form 

u = L-\f - w 2 M 1 L g 1 (L p 1 f) 2 ). (3.12) 
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3.3. Periodic lattice and strip with localised inclusion. In this case the wave equation 
( 13. 7p becomes 

Lu = -w 2 (Mu + Mi(u) 2 + M 2 (u}) (3.13) 
with (■) = ((-)i) 2 . The equation (13.131) is equivalent to 

Lu + w 2 Mu + w 2 Mi(u) 2 = -w 2 M 2 (u). (3.14) 

Applying (13TTTp -( 13~T2|) to (IBTIj) with f = -w 2 M 2 (u) leads to 

u = L- X (i - c^MiL- 1 (L p 1 ) 2 )(-u; 2 M 2 (u)) (3.15) 

and then 

(u) = ^(Lp 1 ),^-^^ 1 ^ 1 )^^-^^)) (3.16) 

which can be rewritten in the form of 

(I + a; 2 (L g 1 (L p 1 ) 2 ) 1 M 2 )(u) = 0. (3.17) 
The equation (I3.17P gives us the condition of existence of localised modes (12. lip . 

4. Example 




a) b) c) 

Figure 2. a) Homogeneous lattice with mass M — 1, b) perturbed by periodic 
strip with mass M = M + M^ l > , c) perturbed by single mass M M • M^ l > + 
M( 2 ). 



4.1. Uniform lattice. We start from homogeneous square lattice with mass M — 1, see Fig. 
|2ja). The wave equation for this structure takes the following form (see (12.31) ) 

(2 cos ki + 2 cos k 2 - A)u = -lo 2 u. (4.1) 
The operator L p and Floquet branches are (see (12. 7p ) 

L p = to 2 - 4 + 2cos/ci - 2cos/c 2 , (4.2) 

w p (k) = a/4 - 2 cos &i - 2cosA; 2 . (4.3) 
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4.2. Uniform lattice with uniform line. Now we add the periodic line with mass 
(recall that M — 1, see Fig. EJb). In this case 



y/(2 cosfci-4+w 2 ) 2 -4' 



v /(2cosfci-4+o; 2 ) 2 -4 

Due to (12. 8p the guided spectrum w g (fci) is determined by 



if oo 2 < 2 — 2cosfci, 
if (jj 2 > 6 — 2 cos k\. 



1 



(det L g 



1 + 



v /(2cosfei-4+aj 2 ) 2 -4 



v /(2cosfci-4+w 2 ) 2 -4 

These equations can be solved directly 



0, if co 2 < 2 - 2cosfci, 
0, if u 2 > 6 - 2cosfci. 



2cosA ;i -4-2 v /(M(i)) 2 (cosfc 1 -2) 2 -(M(i)) 2 +l {1) 

(M(!)) 2 -l ' 11 ^ U ' 

2cosfci-4+2 A /(M( 1 )) 2 (cosfci~2) 2 -(M( 1 )) 2 +l .„ , 

(Mf 1 )) 2 -! ' 11 M > U - 



(4.4) 



(4.5) 



(4.6) 



The projection of the propagative spectrum w p (k) (14.31) on the plane (u, ki) (see (I2.10p ) is 



'l — 2 cos kx, \j 6 — 2 cos fci 



(4.7) 





a) m( 1 )=-o.9 



b) AfW= 



-0.5 



C ) M( 1 )=0.5 



d) tfP)= 



Figure 3. Projection of the Floquet spectrum I p {ki) (14. 7p on the plane (u, k\) 
(red area) and the guided spectrum (green line) cu g (fci) (14. 6p for homogeneous 
square lattice (mass M = 1) with perturbed line (mass see Fie;. 

12b). 



4.3. Uniform lattice and line with localised inclusion. Now we add one perturbed mass 
M = M + M^ 2 \ The condition on uj\ oc (see (12. lip ) takes the following form 



AocM = 1 + 



u 2 M^ 



dki 



2vr J_ n oo 2 MW + ^/(2cosk 1 -A + uj 2 ) 2 -4 



0. 



(4.8) 
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The projection of the propagative and guided spectrum on the axis u (see ( 12.1 3p ) is 

J p = [0,2V2], (4.9) 

2 



V6+2 V /8(M( 1 ))2+1 



0. 



-6+2 v /8(Af< 1 )) 2 +l 
(A/f 1 )) 2 -! 



, if - 1< M« < 0, 
if M W > 0. 



(4.10) 




a) Af( 1 )=-0.9, Af< 2 )=-0.03 b) M< 1 )=-0.9, Af( 2 )=0.1 c) Af( 1 >=-0.9, Af< 2 )=0.25 d) A/< 1 )=-0.9, Af< 2 )=0.7 



e) Af< 1 '=-0.5, A/f( 2 >=-0.2 f) Af< 1 )=-0.5, Af< 2 )=0.1 



I) Af( 1 >=2, M< 2 )=-2.6 



fl) Af«=2, M< 2 )=-2 



Figure 4. Function Di oc (u;) 04.81) (blue curve) and its zeroes u\ oc (localised 
spectrum) with the projection of propagative I p (14. 9 p (red bold line) and guided 
I g (14. 9p (green thin line) spectra on the axe u. 



4.4. Existence of the localised waves (states). The equation ( 14. 8 p can be rewritten in 
the form 

AocH = 1 + M^Dtiuj) = (4.11) 

with 



DM = — / ; 

2tt w 2 M (i) + cos jfcx - 4 + u 2 ) 2 -4 



(4.12) 



The integrand function in (I4.12p is strictly decreasing in a; G K+ \ (J p U J g ) for any k±. Then 
is strictly decreasing function in u G M + \ (I p U J g ). So knowing the values of Di(cu) 
at the edges of I p and/or I g and using (14. lip with monotonicity of Di(uj) we can predict the 
existence of localised modes in the corresponding gap. It is convenient to consider different 
cases (recall that 1 + M^ 1 ' > and 1 + + > because these are the masses of 
inclusions) : 
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4.4.1. The case 6 (—1, — l/y/2). Due to (14 .9p we have two disjoint gaps Gi and G 2 in the 
propagative and guided spectrum, i.e. 

M + \(/ p U/ 6 ) = ftUfe where (4.13) 
( r 2 \ / V /6 + 2 v ^(M( 1 )) 2 + l \ 

g 1 = 2V2, -. g 2 = [ -1 — =^=^ — ,+00 . 4.14 

V \A-(^ (1) )v V Vi-(^ (1) ) 2 / 



The values of Di at the edges of gaps are 

DA2V2) = - / . 1 = < 0, 4.15 

Wo 4M« + ^(008^ + 2)2-1 

/ 2 \ / V / 6 + 2 v ^(M( 1 )) 2 + l\ 
D x — = = -00, £>i( — , =+00, (4.16) 

A(+oc) = TT ^ T y>0. (4.17) 

Using fl4.15p - fl4.17p . (14. lip with monotonicity of D\ we conclude that: 

a) If < then there are no localised states in Gi and there is only one localised state 
in Q 2 . 

b) If = then there are no localised states in Gi and in Q 2 . 

c) If E (0, — l/Z?i(2\/2)) then there is only one localised state in Gi and no localised 
states in Q 2 . 

d) IfM® ^ -l/D x (2y/2) then there are no lacalised states in Gi and Q 2 . 

4.4.2. The case M"W G [— l/y/2, 0]. Due to (14.91) we have only one gap Q 2 in the propagative 
and guided spectrum, i.e. 



/a/6 + 2^8(M( 1 )) 2 + 1 \ 

H + \(j p ui.) = fe = ( v Vl - (M u»p 'H' (418) 

Applying the same arguments as in the previous case we deduce that: 

a) If < then there is only one localised state in Q 2 . 

b) If y> then there are no localised states in Go. 

4.4.3. The case MW > 0. Due to (14.91) we have only one gap G in the propagative and guided 
spectrum, i.e. 

R+\(/ p U/ g ) = £ = (2^+00). (4.19) 

Applying the same arguments as in the previous cases we deduce that: 

a) If < — l/Z?i(2\/2) then there is only one localised state in G- 

b) If ;> -l/ J D 1 (2\/2) then there are no localised states in G ■ 
Note that in this case £>i(2\/2) > (compare with (I4TT5D ). 



4.4.4. Summation. Now we summarize above statements. 
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1) For the strip "I" there are two spectral gaps, 
first gap is located between propagative and 
guided spectra, the second is above the guided 
spectrum (Fig. [3ja). For the brown area the 
localised mode is in the first gap (Figs. HJb,c), 
for the purple one is in the second (Fig. Hja). 

2) For the strip "II" there is one spectral gap 
above the guided spectrum. Guided and prop- 
agative spectra intersect (Fig. Ob). For the 
purple there is one localised mode in the 
gap (Fig. He). 

3) For the quadrant "III" there is one gap 
above the propagative spectrum. Guided 
spectrum lies inside the propagative (Figs. 
Elc,d). For the purple area there is one lo- 
calised mode in the gap (Fig. IUg). 
For the white area we have no localised modes 
in the spectral gaps (Figs. |Ud,f,li). 

Figure 5. The regions of emergence of localised modes for different masses 
M, M in the lattice Fig. EJc). The boundary of these regions consist of curves 
M = M - (see (Q5j) for D x ) and M = M. The maximal mass M for 

which the localised mode exists tends to 1 (with M — > 1). After increasing of 
mass M ^ 1 we see at the beginning the rapid fall of mass M, which reach the 
limit I — it- for M — > +oo. 

5. Conclusion. 

The dispersion equations for the propagative, guided and localised waves in the discrete 
periodic lattice with the strip and with localised inclusions are obtained. For the uniform 
lattice with the line and with one single inclusion the existence of localised modes and its 
position is analysed. 

Acknowledgement. The author is grateful to A. Shuvalov for useful discussions. 
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